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Anomalous enhancement of a penguin hadronic matrix element in B ^ Kr}' 
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We estimate the density matrix element for the 7r° , rj and rj' production from the vacuum in the 
large-A^c hmit. As a consequence, we find that the QCD axial anomaly leads to highly non-trivial 
corrections to the usual flavour SU(3) relations between _B° K°n", B° K°r) and B° K°r)' 
decay amplitudes. These corrections may explain why the B — > Krj' branching ratio is about six 
times larger than the B — > K-r one. 
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First observations of a large branching ratio for B 
Krj' triggered numerous theoretical investigations within 
and beyond the standard model. The latest average of 
the available experimental data Q]: 



Br(B" ^ K^ii) = (64.9 ± 3.5) x 10" 



(1) 



definitely confirms a sizable excess of rj' compared with 
7r°, Br{B^ K°Tr°) = (10.0 ± 0.6) x 10-^ and 77, 
Br{B° K^rj) < 1.9 x 10"^. The fact that rj' is mostly 
a flavour 5'?7(3)-singlet naturally suggests mechanisms 
such as the so-called singlet contribution (i]' production 
by gluon-gluon fusion from b ^ sg or b sgg) or 
intrinsic charm contribution (r/ production by cc annihi- 
lation from b — > ccs) Q. However, these do not seem to 
explain why the B^ — > K^rj' branching ratio is about six 
times larger than the B^ — K^-k^ one. More recently, 
the hadronic parameters for rj' production have been re- 
examined and a couple of different solutions were pro- 
posed HSU. In this letter, we analyse an overlooked 
correction from the axial U{\) anomaly in the hadronic 
matrix elements associated with b sdd, sss. 

The derivative of the flavour-singlet axial current is 
given by: 



q—u,d.s 



47r 



(2) 



where Ga/s is the gluonic field strength tensor and G"'', 
its dual. The mass term in the r.h.s. of Eq. (O implies 
explicit flavour SU{3) violation. This breaking alone 
would lead to the ideal mass relations 



(M2 



'^^ ideal — 2A/|- — Af^, 



(AChdeal = Ml (3) 



which arc quite successful for the ((j),uj,p) vector mesons 
but totally unreahstic for the {r],r]',n) pseudoscalar 
mesons 0. However, the second term in the r.h.s. of Eq. 
(O represents the QCD anomaly which breaks the axial 
[/(I) symmetry to provide the 7/ with a mass around 1 
GeV 
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The hadronic matrix element (OIGq/jG"''!??'^'^) associ- 
ated with the anomalous term in Eq. ^ has been consid- 
ered in relation to the Zweig-supprcssed (J/V' ^ V^'^l) 
radiative decays. In this letter, we aim at an estimate 
of the corresponding hadronic matrix elements for the 
first term in the r.h.s. of Eq. 0. This will then al- 
low us to compute the contribution of the dominant pen- 
guin density-density operator to the {B'-^ K'^n^ ,ri,rj') 
hadronic decays in a way fully consistent with the ax- 
ial U{\) and flavour SU{i) symmetry-breaking require- 
ments. 

In a naive quark picture, the hadronic matrix elements 
(0|g75g|77', 77, 7r°) are simply related through Clebsch- 
Gordan coefficients (C.G.). These coefficients are fixed 
by the qq content in the pseudoscalar wave-functions. 
However, corrections due to flavour SU{?>) but also to 
axial U{1) violations have to be taken into account. In 
particular, we expect the following generic form: 



(0|g759k°) 



C.G.jl 



q—u,( 



A2 



(4) 



The appearance of the physical ry' mass may, at first sight, 
be surprising. But this is in fact required by the ax- 
ial U{1) symmetry. Indeed, in the absence of the axial 
anomaly, the ideal mass relations given in Eq. ^ would 
consistently imply that the rj' /tt^ ratio in Eq Q is equal 
to -fl ii q = u and —1 li q = d since the 7/ and ir^ 
wave-functions have the same quark content as the iso- 
singlct uj and the iso-triplet p in this fictitious world. The 
axial anomaly calls thus for a sizable correction to the 
{Q\q^bq\v') hadronic matrix elements if the cut-off scale 
is what we naturally expect from the real QCD dynamics, 
namely A = 0(1) GeV. So, these matrix elements have 
to be consistently extracted from a low-energy effective 
theory of QCD. 

Again with reference to the observed pattern for the 
pseudoscalar mass spectrum, let us consider the large- 
N(. limit at each order in the (squared) momentum p^. 
The genuine U{Z)l x ?7(3)/f chiral invariant structure of 
the effective non-linear theory implies then the following 
hierarchy Q: 



0{p\ l/N,) > 0{p\ 1/00) > 0{p\ 1/00) . 



(5) 
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such that the full pseudoscalar mass spectrum naturally 
arises in three steps. The leading term ensures 
the breaking of the U{1)a symmetry and provides the 
flavour-singlet 770 with a mass around 1 GeV. The next- 
to-leading (p^) term implies the usual SU (3)y mass split- 
ting among the (tt, K, rjs) flavour-octet. Eventually, next- 
to-ncxt-lcading (p^) terms arc needed to precisely repro- 
duce the observed ?y — rj' mass splitting. 

Starting from this effective theory of QCD in the large- 
Nc limit, we may consistently express the anomalous op- 
erator of Eq. ^ purely in terms of the flavour-singlet 
field rjo. This allows us to extract the rj — rj' mixing angle 
associated with the diagonalization of the rjs — 770 squared 
mass matrix: 



r] = rjg cos 6 — r/Q sin 9 

v' 



?78 sm + rjo cos ( 



(6) 



from the wcU-measurcd {J/tp ^/'■'t) radiative de- 
cays The extracted value is 6* = -(22 ± 1)°. Let 
us emphasise once again that in the absence of the axial 
anomaly, 9 would have been equal to the ideal octct- 
singlet mixing angle -1-35.3° of the <f> — oj, i.e. 



cos 9 ideal 



sm 9 ideal 



(7) 



since the ideal mass relations given in Eq. (jSJ correspond 
to the wave-functions 77 = —ss and if = {uu + dd)/^/2. 
This illustrates how jy'-'-' masses and mixing are strongly 
correlated through the U{1)a symmetry. 

Similarly, we may also express the density operator of 
Eq. (0) in terms of the flavour-nonet field tt: 



-a b ■ f 



1 - 



2V2 



(8) 



A2 



+ r{ma + rrib) 



„ba 



where a,b = 1,2,3 arc the flavour indices. The parame- 
ters r and / are related to the physical masses and decay 
constants as 



M2 = rm,[l + M2(- 



Ml 



and 



/. = /[l + M^(^ + ^)] 
/K = /[l + A/i(^ + ^)] 



(9) 
(10) 



(11) 
(12) 



if isospin symmetry is assumed, q — u 01 d. As a result, 



we obtain the following set of density matrix elements: 
U Ml 



(0|g75(/|7r") = ±i 



2^2 ruq 



. mi^s\K) = ifK 

• (OlgTs?!?/) = 

2^/6 mq 

• (0|g759|?/ 



Ml 



rris + niq 
(c9 - V2s 



A2 
^^0 



2V6 ruq 
(0|s75s|r/) 



{y2c9 



Jk 



K 



Ml - Ml 



M^, - M^ 
A^ 

(V2c9 + s9 



(13) 

(14) 
(15) 

(16) 
(17) 



A2 
^^0 



2(M^-M^)(^-^) 



(18) 



Ml- Ml , , / 1 1 \ 

1+ \, ^ +2{Ml^Ml)[4,--^] 



A2 



A? 4A2; 



where we introduce the abbreviation (c9, s9) for 
(cos0,sin0). On the basis of Eq. we neglect ©(l/A^) 
corrections. This parametrisation of the density matrix 
elements for tt", 77 and r/ follows from a [/(3)l x [/ (3)ij in- 
variant theory and is thus fully consistent with the U{1)a 
and SU{i)v symmetry requirements on the pseudoscalar 
mixing and masses in the isospin limit. 

Eqs. (|13() and H14|l are the well-known hadronic ma- 
trix elements already derived in the large- iV^ limit using 
chiral perturbation theory |lC| . These hadronic matrix 
elements feel the usual effects of SU (3)y violation on the 
decay constants and masses in the pseudoscalar flavour- 
octet. 

Eqs. (|15|l and (|16|) display the highly non-trivial ef- 
fect of the U{1)a breaking: in the absence of the ax- 
ial anomaly, the ideal masses (see Eq. (PJ) and mixing 
(see Eq. Q) would consistently imply (0| 575(7! 77) = 
and {0\q^5q\ri') = ±{Q\q^^q\T:^) . They nicely confirm our 
original guess expressed in Eq. Q). Moreover, a global fit 
of the pseudoscalar masses, mixing and decays constants 
has already been undertaken in our previous work |^. 
The resulting values for the Ao,i,2 cut-offs: 

Ao ~ 1.2 GcV, Ai ~ 1.2 GcV, A2 ~ 1.3 GeV (19) 

are indeed all around one GeV, as anticipated. 

Eqs. (|17|1 and (|18|l consistently combine the effects of 
U{1)a and SU['i)v violations. Our results for the scale- 
independent density matrix elements 2ms{Q\s^zs\'ri^''^) are 
compared with previous works in Table 1. The second 
column includes flavour SU(2i) breaking with a realistic 
octct-singlct mixing angle but with the ideal mass re- 
lations for 77*^'' (sec Eq. (|2l). We can see an excellent 
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this work 


previous works 


e = -22° 


U{1)a X SU{3)v 


SU{3)v 


AG [12J 


BN [4j 


2ims{0 S75S r]) 


+0.053 ±0.008 


±0.058 


±0.057 


±0.055 


2ims(0 S75S 77') 


-0.109 ±0.016 


-0.069 


-0.071 


-0.068 



TABLE I: Comparison of numerical results for the S75S density matrix elements in GeV'^ units. The second column includes 
flavour SU{3) breaking with a realistic octet-singlet mixing angle, 6 = —22°, but with the ideal mass relations for ry''' (see 
Eq. JHJl). Our result agrees with the previous works in this peculiar limit. The first column includes full U{1)a and SU{3)v 
breaking effects. We find that the magnitude of the rj' hadronic matrix element increases by about 60% while the one for rj 
only decreases by about 10%. 



numerical agreement in this limit which is rather pecu- 
liar since any realistic i] — rj' mixing excludes ideal rj — rj' 
masses from the viewpoint of tlie U{1)a symmetry. The 
first column includes full U{1)a and SU{?>)v breaking ef- 
fects. The magnitude of the rj' hadronic matrix element 
increases then by about 60% while the one for rj only de- 
creases by about 10%. Theoretical uncertainties associ- 
ated with higher order SU{^)v and U (1)a corrections are 
expected to be M^/A" ~ 5% and M]^M^^,JA-^ ~ 15%, 
respectively. 

Eqs. H13|) to (|18|) turn out to be crucial for an es- 
timate of the K^rj'-'^ decay amplitudes. Indeed, 
their typical M"^ /m chiral enhancement, at the basis of 
the AI = 1/2 rule in K nn decays is such 
that the Qq = J2q(^LqR){(iRSL) weak operator provides 
the main contribution for the b sqq penguin-induced 
hadronic processes involving two pseudoscalars in the 
final state. (Notice that such is not the case for the 
vector-pseudoscalar final state, B K*rj'''\) In the 
limit where, diagrammatically, B'^ K^rj^'^ come from 
b sdd and b —f sss penguins, while — > K'^tt^ comes 
only from b sdd, we may thus write 



obta 



A{B° K^rj) 
T 



(21) 



-V2 



A{B° K°T]') 

T „ 



-V2 



sex 



(22) 



where the U{1)a and SU{3)v breaking effects associated 
with the pseudoscalar masses are fully encoded in the 
parameters 



c 



(') = 



Ml 



A2 



+ 2(A/ 



1 



1 



^^-^)(a?-4a!)- 



For a sensible value of the mixing angle, 9 = —19.5° 
(i.e. cos 61 = 2^2/3 and sin0 = —1/3), we have rj = 
[uu + dd — ss)/\/3 and rj' ~ (uu + dd + 2ss)/V6, such 
that Eqs. and simply reduce to 



A(B0 ^ ifO^W) 



{K%bLSR\B°] 
{nO\bLdR\BO) 



(20) 



WMbLdRlB") {rj^'^\snSL\0) 
{KO\bLSR\BO) {KO\dnSL\Q) 



A{B° K°Tj) 



A{B^ ^ A'Ovr") 
A{B° K°7j') 



A{BO A'OttO) 



= -19.5° 



?=-19.5o 



-\I-[1-C] (24) 



3[1±2C']- (25) 



We already notice that in the absence of the axial 
anomaly, the ?/ wave-function has no ss component and 
the b — > sss penguin contribution to B^ — > K'^rj' van- 
ishes. In that ideal limit, the rj'/n'^ ratio in Ea. H20() 
is equal to -1 and thus Br{B° -> K^rj') = Br{B" 

As a first order approximation, the effect of the U{1)a 
and SU{3)v violations on the pseudoscalar masses can 
be safely neglected compared to the i3-mass scale. So, 
we simply express the i3-to-light meson hadronic matrix 
elements in terms of C.G. and focus on the axial U{1) and 
fiavour SU (3) breaking corrections for the vacuum-to- 
light meson transitions. In this approximation, we easily 



In the absence of the axial anomaly, the ideal relations 
given in Eq. ® for M^^,) imply (^'^ = 1.41(1.09). If in 
addition we assume Mk ~ M.^, then ( ^ (' = 1 and we 
recover the well-known fiavour SU (3) relations between 
the A'°[0 ''] branching ratios (77' : : vr" = 3 : : 

1) HH' But again, the U{1)a symmetry requires physical 
7/'^ masses for a physical 77 — 77' mixing angle. Hence we 
have 



C = 1.29 ±0.19, C' = 1-72 ±0.26. 



(26) 



The impact of the U{1) anomaly on the — > K'^rj' 
branching ratio is displayed in Fig.l for a mixing angle 
e = -22°. The result with only SU{3)v breaking cor- 
rections (a dot in Fig.l) is obtained for (' — 1.09. As we 
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FIG. 1: The B" K^-q' to B" K°tt° ratio of branching ratios as a function of the 77' mass-dependent parameter (' , for a 
fixed value of the mixing angle 6 = -22°. The usual flavour SU{3) relation, Br{B° K°r]')/Br{B° K°n") ~ 3, is obtained 
at C — 1. The dot corresponds to SU{3)v but no U{1)a corrections to (' , i.e. — 1.09. Finally, the star represents full 
517(3)^ and U{1)a corrections to (' , i.e. = 1.72 ± 0.26. The shaded area displays the current experimental range. 



have shown in Tabic 1 , this pecuUar hmit reproduces the 
numerical values of the hadronic matrix elements used 
in the previous works. Our result which consistently in- 
cludes the effect of the U {1)a violation on the masses and 
mixing (a star in Fig.l) is based on Eq. H26() for C'. We 
observe a strong increase for the 5° K'^rj' branching 
ratio, in agreement with the data. On the other hand, the 
relatively small 77 mass induces a small U{1)a correction 
to which does deviate from one mainly through the 
SU{3)v violation. As a result, the B'^ K'^rj branch- 
ing ratio stays well below the experimental bound due to 
the efficient cancellation in Eq. H24(l but is, at the same 
time, extremely sensitive to the uncertainties on ^ and 
e. For C = 1-29 ± 0.19 and = -(22 ± 1)°, its value 
runs between 0.01 x 10"*^ and 1.10 x 10"^. At this level, 
contributions from the other penguin operators might be 
non-negligible. 

In conclusion, we reconsidered the puzzle of large 
B Krj' branching ratio and pointed out a missing 
U{1)a breaking correction to the penguin hadronic ma- 
trix element. We estimated this correction in the frame- 



work of a low-energy effective theory of QCD in the large- 
Nc limit. We provided the expression for all the density 
matrix elements relevant to the hadronic B decays, in 
terms of physical masses and decay constants. We found 
a rather large increase (60%) of |(0|s75s|77') | and a mod- 
erate decrease (10%) of |(0|s75s|7y)| compared to the pre- 
vious works. The sizable corrections for the rj' should 
not come as a surprise since its hadronic matrix element 
vanishes in the absence of the axial anomaly. This cor- 
rection may explain why the B'^ K'^rj' branching ratio 
is about six times larger than the _B" K'^tt'^ one. 
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